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Abstract.  The  validity  of  a  common  radiometeorological  application  of  Monin-Obukhov  (M-0)  simi¬ 
larity  theory  to  potential  refractivity  (;y),  which  is  a  nonlinear  combination  of  d  and  q,  is  determined 
by  whether  the  properly  nondimensionalized  gradient  is  a  universal  function  of  zIL.  We  develop 
expressions  for  the  flux  of  y  (and  its  "scaling  parameter,  x*)  terms  of  temperature  and  moisture 
fluxes,  and  an  M-0  similarity  expression  for  the  vertieal  x  gradient.  Results  show  that  even  if  9  and 
q  are  accepted  as  exactly  following  M-0  similarity  expressions;  when  the  surface  layer  is  stable,  x 
does  not  obey  such  an  expression.  That  is,  when  properly  nondimensionalized,  the  vertical  gradient 
of  X  does  not  collapse  to  a  single  universal  function  of  zIL,  i  he  assumption  that  x  behaves  as  a 
similarity  variable  is  approximately  correct  for  well-mixed  surface  layers  under  unstable  and  near¬ 
neutral  conditions. 

The  gradient  of  x  is  an  important  factor  in  determining  microwave  propagation  conditions.  We 
demonstrate  the  error  induced  in  a  simple  algorithm  when  x  is  assumed  to  obey  M-0  similarity  theory. 
An  alternative  methodology,  consistent  with  the  application  of  similarity  theory  to  9  and  q,  is  then 
developed  without  requiring  that  x  itself  satisfy  similarity  theory. 


1.  Introduction 


Microwave  propagation  near  the  sea  surface  depends  primarily  on  the  refractive 
structure  of  the  marine  atmospheric  boundary  layer.  The  refractivity,  or  refractive 
index,  of  air  is  related  to  its  density  and  is  a  nonlinear  function  of  temperature 
and  humidity.  Analysis  and  modeling  of  the  mean  temperature  and  moisture 
structures  within  the  surface  layer  have  a  rich  history,  resulting  in  what  are 
generally  referred  to  as  flux-profile  relationships  based  on  Monin-Obukhov  (M- 
O)  similarity  theory  (Monin  and  Obukhov,  1954).  Since  refractivity  is  so  closely 
related  to  these  similarity  variables,  it  has  been  widely  assumed  that  it  can  be 
treated  as  a  similarity  variable  as  well  (e.g.,  Anderson  and  Gossard,  1953;  Gossard, 
1964;  Brocks,  1965;  Jeske,  1971;  Hitney,  1975;  Gossard,  1978;  Gossard,  1981; 
Thompson,  1987;  and  Paulus,  1989).  In  particular,  these  authors  all  use  similarity 
relationships  involving  potential  refractivity  (;^): 


APo 
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1  + 


Bq 

€d 


(1) 


where,  for  microwaves  with  wavelengths  greater  than  1  cm,  A  and  B  are  constants, 
equal  to  77.6 °K mb"*  and  4810 ®K,  respectively  (Bean  and  Dutton.  1966),  Po  = 
1000  mb,  6  is  the  potential  temperature,  q  is  the  specific  humidity,  and  e  =  0.622. 
Potential  refractivity,  rather  than  normal  refractivity,  is  generally  used  because 
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it  has  been  felt  that.as  a  conservative* scalar  (like  d  and  <y);  it  necessarily  must  also 
be  a  gcod  variable  for  similarity  theory  applications.  We  shall  demonstrate  how¬ 
ever,  that  ;k' cannot  be  treated  as  a  true  similarity  variable  in  a  scheme  consistent 
with  the  treatment  of  B  and  q.  We  shall  assess  the  magnitude  of  the  error  one 
makes  when  assuming  that  similarity  theory  applies  to  x  atid  show  that,  for 
practical  purposes,  x  can  be  tfeated'as  a  similarity  variable  with  little. error  under 
unstable  and  near-neutral  conditions.  This  may  explain  its  widespread  usage  in 
this  manner.  Under  stable  conditions,  however,  x  deviates  sufficiently  from  simi¬ 
larity  theory  to  cause  significant  errors  in  its  application. 


2.  Methodology 

Using  Reynolds  averaging,  the  turbulent  flux  of  a  scalar  quantity,/,  at  the  surface 
is  (Fleagle  and  Businger,  1980): 

where  p  is  the  surface  layer  atmospheric  density,  w'  is  the  turbulent  component 
of  the  vertical  wind,  and  /'  is  the  turbulent  fluctuation  of  the  scalar  /.  M-0 
similarity  theory  states  that  in  the  surface  layer,  the  vertical  gradient  of  such  a 
scalar  quantity  is  independent  of  the  character  of  the  surface  (e.g.,  roughness) 
and  is  determined  by  the  density  p,  the  “buoyancy  parameter”  glBa,  a  measure 
of  the  turbulent  shear  stress  u*,  a  measure  of  the  vertical  heat  flux  d*,  and  £/, 
the  flux  of  /  (see  Monin  and  Yaglom,  1971,  Section  7.2,  for  a  discussion  of 
why  these  particular  parameters  are  relevant).  Furthermore,  the  gradient  may  be 
described  by 


1^1 
/()  Bz 


(2) 


where  F(zlL)  is  a  universal  function  of  z/L.  The  Obukhov  length  is  defined  by 
L  =  iddalgkd*  where  g  is  the  acceleration  due  to  gravity  and  k  is  von  Kerman’s 
constant.  Here/o  has  the  dimensions  of  df/dz  and  based  on  the  dimensions  of  the 
available  parameters,  we  may  write  (Monin  and  Yaglom,  1971): 


/o  = 


pu*L  ’ 


and  Equation  (2)  becomes 


L  y 

f*  Bz 


=  FU) . 


(3a) 


where  C  =  zIL.  The  scaling  parameter, /sc ,  is  defined  by 
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In  its  more  standard  form,  Equation  (3a)  is  usually  written 


kaz  df_  __  kaz 
f:lf  dZ  Li 


F(o  ^  m , 


(3b) 


(4) 


where  «I>(^)  is  another  universal  function  and  a  is  the  frtio  of  the  turbulent 
diffusivities  of  the  scalar  property  /  to  momentum  at  neutral  stability. 

In  the  surface  layer,  nondithensional  expressions  for  the  potential  temperature 
and  specific  humidity  gradients  have  been  derived  from  siniilarity  theory  (Moniii 
and  Obukhov,  1954;  Businger  et  al.,  1971),  and  may  be  written 


kazctQ 

{5a) 

Dz 

(5b) 

where  0*  and  q*  are  scaling  parameters  and  0(^)  is  the  universal  function  defined 

by 


<t>U)=l  +  bj-,  (5c) 

for  stable  conditions  where  =  7,  and 

-1/2 

m  =  ,  (Sd) 


for  unstable  conditions  where  a  =  16  (Liu  et  al.,  1979). 

When  the  scalar  property  /  equals  potential  refractivity  we  wish  to  know 
whether  Equation  (4)  is  satisfied.  However,  to  test  Equation  (4)  we  need  the 
surface  potential  refractivity  flux,  E^.  We  now  describe  a  method  for  relating  the 
potential  refractivity  flux  to  the  temperature  and  moisture  fluxes. 

If  we  assume  that  Equation  (1)  applies  to  instantaneous  values  of  and  we 
decompose  the  dependent  variables  into  means  and  turbulent  perturbations  (e.g., 
A"  =  ^  +  A''>  etc.),  then  we  may  write 


Expanding  binomially,  neglecting  second-order  and  higher-order  terms,  and  sub¬ 
tracting  the  mean,  x,  yields 


>0 
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Multiplying  by  pw'  and  performing  a  Reynolds  average  gives' the-flux  of  as 
r-  .—7^^  /4Po  .  24flPo^^  ^-=7X7  /4BPo\ 


and  thus,  from  Equation  (3b),  the  scaling  parameter  to  be  iisbd  in  Equation  (4) 
is 


Ar*  = 


\  e  €6  J  \  €d  - 


q* . 


(6) 


We  note  that  this  scaling  parameter  can  also  be  derived  by  taking  the  vertical 
derivative  of  Equation  (1)  and  substituting  Equation  (5)  for  Bdlbz  and  dqidz. 

To  determine  whether  x  behaves  as  a  similarity  variable,  one  would  ideally 
use  observational  data  to  verify  Equation  (4).  Unfortunately,  accurate,  abundant 
marine  observations  of  x  are  not  available  so  we  developed  a  numerical  procedure 
to  test  the  relationship.  This  approach  has  the  advantage  of  eliminating  consider¬ 
ation  of  random  observational  errors  from  the  analysis.  In  our  procedure,  the 
integrated  form  of  Equation  (5)  from  Liu  et  al.  (1979)  was  used  to  compute 
similarity  profiles  of  B  and  q  at  Az  =0.im  intervals  from  0  to  40  m  above  the 
surface,  from  which  we  obtain  x  using  Equation  (1).  BxiBz  at  each  level  was  then 
calculated  using  centered  finite  differences  with  truncation  error  of  0(Az‘’).  We 
calculate  x*  by  evaluating  Equation  (6)  using  B  and  q  at  the  reference  height, 
taken  in  this  paper  to  be  2  m.  Then  we  compute  the  left-hand-side  of  Equation 
(4),  with  f=  x^  at  each  grid  point  for  a  wide  range  of  stability.  The  results  are 
then  plotted  as  a  function  of  ^  to  see  if  the  nondimensionalized  x  gradient  is  a 
universal  function  of 


3.  Results 

Figure  1  shows  the  characteristic  surface-layer  profile  curvatures  of  A0,  At?,  and 
A;^,  which  were  computed  from  Equations  (5)  and  (1)  for  both  stable  and  unstable 
conditions.  Difference  quantities  were  chosen  to  eliminate  scale  differences  be¬ 
tween  the  stable  and  unstable  cases.  The  unstable  profiles  in  Figure  la  were 
computed  with  arbitrarily  chosen  values  of  sea  surface  temperature  equal  to  18  °C, 
and  reference  height  values  of  air  temperature  equal  to  12  °C,  wind  speed  equal 
to  3.6ms~‘,  and  relative  humidity  equal  to  80%.  The  pressure  at  the  reference 
height  was  set  at  1000  mb  and  was  computed  hydrostatically  away  from  the  refer¬ 
ence  height.  These  profiles  show  characteristic  logarithmic  shapes  reflecting  sharp 
gradients  near  the  sea  surface  and  nearly  constant  values  up  into  the  well-mixed 
surface  layer.  Under  stable  conditions,  such  as  those  found  near  the  Gulf  Stream 
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UNSTABLE  STABLE 

-  -  -  Specific  Humidity  Difference  (g/kg) 


Potential  Temperature  Difference 


I  I  I  I  I  I _ I _ I _ I _ I _ i _ I _ 1 — I — 1 — I — I 

-48  -42  -36  -30  -24  -18  -12  -6  0  6  12  18  24  30  36  42  48 

-  Potential  Refractivity  Difference 

Fig.  1 .  (a)  Profiles  under  unstable  conditions  of  potential  temperature,  specific  humidity,  and  potential 
refractivity  difference  between  the  profile  at  height  z  and  the  sea  surface.  Computed  for  sea  surface 
temperature  equal  to  18  °C,  and  reference  height  values  of  air  temperature  equal  to  12  "C,  relative 
humidity  equal  to  80%,  and  wind  speed  equal  to  3.6m  s~*.  (b)  Same  as  Figure  la  except  for  stable 
conditions  with  air  temperature  equal  to  24  °C. 

(Ross,  1981;  Hayes,  1981)  or  in  the  North  Sea  (Keller  et  al.,  1989)  when  warm 
air  advects  over  cooler  water,  the  6  and  q  profiles  in  Figure  lb  show  surface-based 
inversions  with  constant  slopes  above  a  shallow  interfacial  layer;  the  slope  of  the 
X  profile  is,  however,  not  constant.  These  profiles  were  computed  with  the  same 
parameters  as  above  except  for  reference  height  air  temperature  equal  to  24  °C. 
Since  a:  is  a  nonlinear  function  of  6  and  q,  its  vertical  gradient  is  not  constant  even 
though  the  vertical  gradients  of  $  and  q  are.  Under  unstable  conditions,  6  and  q 
are  nearly  constant  with  height,  and  thus  so  is  X‘ 

Figure  2  shows  <I>;^(^)  as  computed  from  Equation  (4)  with  /=  x->  and  <^(^) 
from  Equation  (5)  for  many  different  profiles  over  a  wide  range  of  stability.  The 
reference-height  air-sea  temperature  difference  was  varied  from  -1.5°  to  +6°C 
in  0.25  °C  steps  white  sea  surface  temperature,  wind  speed,  and  relative  humidity 
wei“  fixed  at  the  values  cited  above.  This  figure  shows  that  under  stable  conditions. 
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Fig.  2.  Profiles  of  the  <P;^(f)  (squares)  and  <!>(()  (solid  line)  for  near-neutral  and  stable  conditions. 
The  environmental  conditions  correspond  with  those  in  Figure  1  except  here  the  air-sea  temperature 
difference  has  been  varied  from  - 1.5  °  to  +6'C  in  0.25  *C  increments.  The  lines  made  by  superposition 
of  the  squares  are  members  of  a  family  of  curves. 


<I>;f  is  clearly  not  a  universal  function  of  that  is,  the  curves  do  not  fall  on 
top  of  one  another.  <1>;^  and  (ft  are  nearly  coincident  under  unstable  and  near 
neutral  conditions. 

Although  is  a  function  of  two  variables,  6  and  q,  which  satisfy  M-0  similarity 
theory,  x  itself  does  not  satisfy  the  M-O  hypothesis  [Equation  (4)]  because  of  its 
nonlinear  functional  dependence.  Virtual  potential  temperature,  6,.,  is  another 
nonlinear  combination  of  0  and  q  to  which  similarity  theory  is  frequently  applied. 
To  investigate  whether  this  assumption  is  justified,  we  repeated  our  experiment, 
computing  0,.  instead  of  x-  Our  results  (not  shown)  verify  that  0,.  closely  follows 
a  universal  profile  under  all  of  the  environmental  conditions  we  investigated. 

Applications  of  as  a  similarity  variable  are  often  used  to  assess  the  height  (S) 
of  a  critical  potential  refractivity  gradient  (dxldz)c  by  using  Equation  (4)  with  z 
replaced  with  8,  dx/dz  replaced  with  (dx/dz)c,  and  assuming  (Anderson 

and  Gossard,  1953;  Gossard,  1964;  Brocks,  1965;  Jeske,  1971;  Hitney,  1975; 
Gossard,  1978;  Gossard,  1981;  Thompson,  1987;  and  Paulus,  1989): 
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Fig.  3.  Height,  5,  of  the  critical  gradient  (dxldz)^  versus  stability  (f),  calculated  using  the  direct 
method  (solid  line)  and  assuming  x  follows  M-0  similarity  theory  (dashed  line).  The  two  different 
techniques  are  described  in  the  text.  Computed  values  of  5  greater  than  40  m  have  been  set  to  40  m 
and  values  of  S  less  than  Om  have  been  set  to  Om. 


kaS  / _  ,/j5\ 
\L/ 


(7) 


The  height  8,  usually  called  the  evaporation  duct  height,  is  an  important  factor  in 
determining  shipboard  radio  and  radar  propagation  conditions.  Generally,  x*  is 
evaluated  from  the  integrated  form  of  Equation  (4)  using  bulk  parameter  values 
at  a  reference  height,  and  then  Equation  (7)  is  solved  for  5,  Figure  3  is  a  plot  of 
8  versus  stability  (^)  for  a  subset  of  the  cases  presented  in  Figure  2.  The  dashed 
line  represents  the  solution  to  Equation  (7)  and  the  solid  line  (labeled  ‘direct 
method’)  represents  the  solution  to  the  integrated  form  of  Equation  (5),  with  the 
computation  of  9,  q,  x*  and  dxidz  directly  at  each  profile  level,  and  a  searching 
algorithm  to  find  the  height  where  the  x  gradient  reaches  its  critical  value.  Thus, 
the  direct  method  only  assumes  that  6  and  q  satisfy  M-0  similarity  theory  -  not 
that  X  does.  In  Figure  3,  the  direct  method  shows  that  nowhere  in  the  profile, 
from  about  0.1  <^<0.8,  does  the  x  gradient  reach  its  critical  value.  (In  the 
figure,  the  value  of  8  has  been  set  to  40  m  when  calculated  values  exceed  40  m, 
and  set  to  zero  when  negative  values  for  5  are  calculated.)  However,  the  solutions 
to  Equation  (7),  which  do  assume  that  x  follows  similarity  theory,  give  erroneously 


>v, 
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low  values  for  S  in  these  cases.  Other  cases  (hot  shown)  show  that  application  bf 
Equation  (7)  oh  the  stable  side  lihder  different  ehvifohmehtal  conditions  cari  also 
lead  to  erroneously  large  values  of  5;  thus  eveh  the  sign  of  the  effor  that  one 
makes  wheh  treating  x  in  this  riianhef  is  hot  consistent;  This  sign,  change  cah  also 
be  inferred  from  Figure  2  because  can  be  either  smallef  of  larger  than  0  for 
the  same  stability. 

In  addition^  we  found  the  magnitude  of  {dxldi)^  to  have  up  to  10%  vafiaiioh 
with  height  (not  shown);  This  variation  vvas  taken  into  account  in  the  direct  method 
but,  for  reasons  of  Computational  simplicity,  a  constant  value  of  (dx/dz)c  = 
-  0.131  m“‘  was  used  when  solving  Equation  (7).  The  use  of  a  constant  value  for 
(dxldz)c  in  Equation  (7)  is  consistent  with  operational  Navy  computer  codes  (e.g;, 
Paulus,  1989)  and  is  a  small  additional  source  of  error  in  the  estimation  of  5. 

4.  Concluding  Remarks 

We  have  shown  that  even  if  6  and  q  are  accepted  as  exactly  following  M-0 
similarity  expressions,  when  the  surface  layer  is  stable,  x>  a  nonlinear  combination 
of  d  and  q,  does  not  obey  such  an  expression.  That  is,  when  properly  non- 
dimensionalized,  the  x  gradient  does  not  collapse  to  a  single  universal  function  of 
zIL.  For  practical  purposes,  we  found  another  nonlinear  combination  of  B  and 
q,  to  obey  similarity  theory.  It  is  the  functional  form  of  x  that  is  important  in 
determining  whether  or  not  it  obeys  similarity  theory.  When  the  surface  layer  is 
unstable  or  near-neutral,  x  closely  follows  similarity  theory.  It  was,  therefore, 
necessary  to  perform  the  full  range  of  numerical  calculations  that  we  present  here, 
rather  than  try  to  make  a  judgement  about  as  a  M-0  variable  based  solely  on 
its  functional  form.  The  advantage  of  our  numerical  methodology  is  that  the  results 
in  Figure  2  can  be  interpreted  as  systematic  errors  due  to  inappropriate  application 
of  M-0  similarity  theory  to  x-  This  source  of  error  needs  to  be  considered  in 
addition  to  any  error  due  to  measurement  scatter.  Ordinary  ref-activity  is  even 
less  suitable  as  a  similarity  variable  due  to  its  additional  functional  dependence 
on  pressure. 

The  methodology  and  algorithms  embodied  in  Hitney  (1975)  and  Paulus  (1989), 
which  assume  that  obeys  similarity  theory,  have  been  used  in  military  operational 
computer  codes  for  over  ten  years.  Results  from  these  codes  are  known  to  be 
erroneous  under  certain  circumstances,  but  the  corrections  (Paulus,  1985)  have 
not  addressed  the  problem  we  discuss  here.  The  straightforward  alternative  to 
treating  x  as  an  M-0  similarity  variable  is  to  apply  similarity  theory  to  B  and  q, 
and  then  to  compute  refractivity  as  v?  do  in  our  direct  method. 
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